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ABSTRACT 


The air-bubbling technique utilized by the Military Sea 
Transportation Service to prevent ice formation in North Star 
Bugt, Thule, Greenland is discussed. A proposed model of 
induced water circulation is presented. Physical processes 
impeding ice formation and growth in sea water are described. 
Oceanographic data collected in conjunction with the MSTS 
project are analyzed and presented in the appendixes. 


FOREWORD 


The formation of sea ice in northern waters often hastens the 
termination of shipping at Arctic harbors and sites. Retardation of 
ice growth can reduce or prevent damage to marine structures such as 
De Long pier at Thule, Greenland. A thorough understanding of physi- 
eal effects which delay ice formation and slow ice growth is required. 
For these reasons, the experiment performed at Thule in 1959 is of 
considerable interest. 


This report is a study of the effects of air-bubbling on the physi- 
eal properties of the water adjacent to De Long pier. It attempts to 
formulate a working hypothesis for explaining the mechanism of the proc- 
esses which retarded formation and growth of sea ice. 


Conclusions expressed in this report may require revisions as 
additional data become available. All additional information which 
might amplify or modify this report will be welcomed by the Hydro- 
graphic Office. 
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INTRODUCTION 


Successful application of an air bubbling technique for prevent- 
ing ice formation during the fall of 1958 enabled the Military Sea 
Transportation Service to extend the shipping season at Thule, Green- 
land. The normal shipping season extends from early July to early 
October. Shipping during the first half of July is usually dependent 
on icebreaker escort; shipping is ordinarily terminated prior to 
initial ice formation in autumn, 


Adaptation of a method originally developed in Scandinavian coun- 
tries for prevention of fresh-water freezing permitted maintenance of 
an ice-free area (polynya) adjacent to De Long pier despite normal 
ice formation in the surrounding waters of North Star Bugt (Bay). 
Safeguarded against becoming frozen-in at the pier, ships of the 
supply convoy remained at Thule until 25 October 1958 - the latest 
date on which MSTS had ever operated in such a northerly location. 
Success of the temporary installation prompted the Commander Military 
Sea Transportation Service Atlantic (COMSTSLANT) to formulate plans 
for the establishment of a permanently installed bubbling system at 
De Long pier. 


In the fall of 1959, the U. S. Navy Hydrographic Office was 
requested to conduct oceanographic studies concurrently with the 
operation of the system in order to obtain information on physical 
processes impeding ice formation and growth in the bay. The over- 
all operation was dubbed "Project Polynya". 


NORTH STAR BUGT 


North Star Bugt, approximately 3 square miles in area, recedes 
about 1-1/2 miles northeastward between Astro Pynt and Mount Dundas 
on the southern shore of Wolstenholme Fjord. The entrance of the 
bay, about 3/4 mile wide, is narrowed by De Long pier and a causeway 
which extend approximately 0.4 mile west-northwestward from Astro 
Pynt. The pier, 1,000 feet in length and 50 feet in width, adjoins 
the causeway and is parallel to it. Inside the entrance, the width 
of the bay increases to about 1 mile. 


POLYNYA INSTALLATION 1959 


Under the supervision of MSTSLANT, Canadian Underwater Demolition 
Unit BRAVO began installation of the bubbling system in North Star Bugt 
early in September 1959. Briefly, the system consists of perforated 
submarine polyethylene pipes which serve to conduct compressed air to 
the bottom of the bay and distribute it in the form of bubbles over a 
wide area (see Figure 1). An auxiliary alcohol-injection system is 
utilized to prevent or eliminate ice formation within the pipes. 


Installation was completed on 27 September. Equipment was inter- 
mittently tested until appreciable ice began to form in the bay on 
T October, at which time full operation was begun. 


ICE FORMATION AND GROWTH IN 1959 


Ice was first observed on 16 September along the eastern fringe 
of North Star Bugt at a point where fresh water flows from the Pitufik 
River. By 2h September grease ice formed in an area northeast of the 
pier, where shallow water was observed to cool to the freezing point 
during low tide. This grease ice drifted westward with the next ebb 
tide. By 7 October a considerable amount of pancake ice had formed in 
the area northeast of the pier. 


During the morning of 8 October a foehn wind arose with gusts to 
approximately 50 knots. The air temperature increased from 17° to 
340 F. On the following morning, the bay was completely clear of ice 
and the surface water temperature had increased from-1.22° to -0.949 ca. 


During the morning of 10 October grease and pancake ice of small 
‘diameter formed over most of the area northeast of the plier. On 11 
October a considerable amount of large pancake ice was observed drift- 
ing into the region from the south. The areas to the south, west, and 
northwest of the pier attained a coverage of approximately nine-tenths3; 
no major ice formation was observed near the pier. By 15 October most 
of the bay was covered with young floes averaging approximately 10 feet 
in diameter. Yen-tenths concentration of young ice accumulated every- 
where by 20 October except in the ice-free area adjacent to the pier. 


DATA COLLECTION 


Serial temperature, salinity, and current data were obtained at 
the locations shown in Figure 2. Temperature and salinity observa- 
tions were taken between 9 September and 21 October. 


Station 1 was occupied daily using a small hand winch mounted on 
the pier or by use of an oceanographic winch aboard the WESTWIND. For 
detailed observations, the polynya was divided into Stations 1A, 1B, 
and 1C as shown in Figure 3. Station 1D was designated at a point ly- 
ing approximately 15 feet outside the polynya area and near the east- 
ern end of the pier. Stations 2, 3, and } were occupied twice weekly 
using either a Greenland cruiser or an LCVP. Station 5 was occupied 
only on 15 September and Station 1D only on 13 October. Stations 6, 
7, 8, and 9 were occupied weekly by the WESTWIND. 


A 200-foot bathythermograph was used to measure temperature at 
all oceanographic stations. BT drops were made daily within the 
polynya area. Original plans had included daily observations at Sta- 
tions 1, 1A, and 1B; however, presence of shipping often precluded 
data-collection at all 3 locations. A bucket thermometer was used in 
conjunction with each PT drop to obtain surface water temperature. 


Meteorological observations obtained daily on the pier include 
wind speed and direction, wet and dry bulb air temperature, cloud 
cover and type, visibility, sea state, and state of weather. In 
addition, pertinent data were extracted from the weather log main- 
tained at Thule Air Base for analysis (Appendix V). Data on tidal 
currents were obtained by suspending 3 Roberts radio current meters 
at depths of 6, 26, and 54 feet from an anchored radio buoy. Water 
depth at each location was 12 fathoms. Signals transmitted from 
the buoy were recorded at 30-minute intervals by a monitoring sta- 
tion aboard the WESTWIND. Current meter Station 10 was occupied from 
1830Z, 25 September to 2100Z, 26 September; Station 11 was occupied 
from 1200Z, 6 October to 1930Z, 13 October 1959. Current data were 
not tabulated, because the recorded results either approximated the 
threshold value of the current meter (0.2 knot) or were unreadable. 


A supplementary survey conducted during April 1960 by Hydro- 
graphic Office ice observers yielded late-winter temperature data 
at Stations 1 and 2 (Appendix III). 


DATA ANALYSIS 


Oceanographic conditions at each station were examined for 
factors contributing to the formation and growth of ice. Surface 
temperatures indicating heat loss at the sea surface and physical 
properties showing the distribution of heat loss throughout the 
water column were studied. Data obtained outside the polynya 
were compared to those obtained at Station 1 in order to deter- 
mine the effect of the bubbling system on the oceanographic struc- 
ture. 


Reversal of the heat budget had occurred prior to inception 
of the oceanographic survey. Except for interruption by the foehn 
on 8 October, progressive cooling was observed at all depths. The 
temperature rise shown by the tabulated data for 9, 10, and 11 Octo- 
ber was observed at Station 1. Upon resumption of the cooling pro- 
cess, surface temperature outside the bubbled area decreased rapidly. 
The freezing point was attained on 13 October. 


A study of the salinity structure indicates spatial and tempo- 
ral fluctuations of surface values and depth of the isohaline layer. 
A plot of the surface salinity values at Station 1 is presented in 
Figure 4. The portion of the plot constructed from values for early 
September indicates that a certain periodicity may exist. The pro- 
nounced increase during the latter part of September is attributed 
to cessation of runoff. 


Data obtained at Station 3 on 18 September and 12 October are 
plotted in Figures 5 and 6. The surface water temperature on 18 
September was 0.649 Cs; temperature maximum of 0.72° C occurred at 
9 and 20 meters. The surface salinity was 31.23 O/o0. Convection 
extended to a depth of only a few meters. 


By 12 October the surface water had cooled to -1.32° C3; the warm- 
est water was at the bottom. ‘The temperature maximum of -0.62° C was 
observed at a depth of 33 meters. The surface salinity had increased 
to 32.34 °/oo, and convective mixing had produced an isohaline layer 
in the upper 15 meters. The calculated freezing point of the surface 
as EILGyOl Oa 


The density gradient below the 15-meter level, although weak, has 
significant relevance to the bubbling system. Theoretical ice-poten- 
tial calculations using the data of the deeper stations show that, 
prior to initial ice formation , thermohaline convection takes place 
to a depth of approximately 15 meters. Consequently, the water below 
this level temporarily serves as a source of sensible heat. However, 
as ice forms, the salinity of the upper layer increases, resulting in 
greater density and an increase in depth of convective mixing. As the 
density gradient weakens and eventually disappears, cooling to the 
freezing point will occur throughout the water column. At the known 
average rate of heat loss from the sea surface in the latitude of Thule, 
the entire supply of warm water will be eliminated within two weeks 
after initial formation of ice. 


The heat content of air issuing from the compressors is consider- 
able. An appreciable amount of heat.is possibly introduced into the 
bottom water adjacent to the pier when compressed air cools in the 
polyethylene pipes; however, the data do not indicate a temperature 
differential attributable to this source within the bubbler field. A 
layer of dirt covering four steel feed pipes provides insulation; how- 
ever, heat loss through the rubber feed hoses is great. Between the 
point where these hoses connect to the steel pipes and the point where 
they enter the water, melting of snow within a radius of approximately 
2 feet was observed. Ice formation due to moisture condensate in the 
underwater sections of the feed pipes was removed by alcohol injection. 


The effectiveness of the bubbler system, when warmer bottom water 
is available, is manifested by the temperature data in Table I. 


TABLE I 
13 October 1959 

STATION 1 STATION 1D 
Depth Tempe Depth Temp e 
(meters) (9) (meters ) (2c) 

0.0 ~1.09 0 -1.76 

365 =OG 3 -1.60 

6.5 -1.06 6 -1.25 

965 -1.06 9 -1.24 


The surface temperature at Station 1D, located immediately outside 
the bubbler field, shows that the surface water had cooled to the 
freezing point. Despite ice formation around the perimeter of the agita- 
ted area, the data obtained at Station 1 show the surface temperature to 
be 0.67° C above the freezing point. ‘The data for Station 1D are assumed 
to be indicative of temperature data that would have been observed at 
Station 1 had the bubbler system not been in operation. ‘The data of 13 
October plus the profiles for Stations 2 and 3 on 12 October indicate 
that water from depths greater than 15 meters is circulated into the 
agitated water columns adjacent to the pier. 


Proof that the bubble system acts as a huge pump capable of perform- 
ing work on the surrounding subsurface water is provided by comparison 
of data presented in Tables II and III. 


TABLE II 
11 October 1959 
STATION 1 
Depth Temp Salinity e 
(meters) (2G) (°/d0) t 
O50 -0.90 32,42 26.08 
B55) -0.91L 32.42 26.08 
6.5 -0.91 32.42 26.08 
9.5 -0.89 32.43 26.09 


TABLE III 
12 October 1959 
STATION 2 
Depvch Temp Salinity 
(meters) (2 c) (2/00) L 
@) =i .42 32,32 26.02 
5 -1.40 32632 26.02 
15 = .42 32.632 26.02 
22 -0.69 32.52 26.16 
eh -0.78 32.52 26-16 


The density of the agitated water column in Table II is greater 
than the density to at least 15 meters in Table III; therefore, work 
was performed by the system in raising water through a vertical dis- 
tance in excess of 15 meters. Comparison of salinity and density 
data of Table II and the plotted curves of Figure 6 reveals that 
water similar to the entire water column at Station 1 is found at 20 
meters at Station 3, indicating that the water was raised at least 
20 meters. 


The eventual cooling of the entire water column to the freezing 
point indicates that vertical transport of sensible heat from depth 
was not a factor in the maintenance of the artificially created 
polynya, except during the initial stage of the ice formation. Con- 
sequently, an understanding of the physical process involved must be 
sought along other lines. 


Elementary ice particles are probably disk-shaped and devoid of 
erystalline form. Ordinarily they flocculate and grow into true 
erystals. The turbulent energy of the induced currents may destroy 
the crystals before they enlarge or may effectively prevent crystal- 
line growth about ice nuclei. Ice particles at the surface of the 
bubbled area are rapidly swept from regions of divergence into 
regions of convergence where, by means of descending currents, they 
are transported beneath the surface to be eventually dispersed from 
the polynya area. 


Hydrographic Office ice observers, stationed at Thule Air Base 
throughout the winter of 1959-60, noted that the polynya gradually 
narrowed; by the end of December width ranged from 12 feet at the 
eastern end to 50 feet at the western end, where an auxiliary air 
hose was used to augment the bubbling activity by inducing more vig- 
orous currents. Dimensions of the ice-free area gradually increased 
during spring as the air temperature rose to approximately 0° F. 


A plot of sea ice tensile strength versus temperature (Assur, 
1958) shows a marked increase of strength as the temperature of the 
ice drops below -9.2° F. At this temperature sodium chloride is 
precipitated from the brine pockets in the ice. During periods of 
extremely low air temperature in winter, the weakest point of the 
ice should be at its undersurface where the temperature approaches 
that of the water. 


Measurements made during April 1960 show that ice thickness 
directly above one of the polyethylene pipes averaged approximately 
10 inches while thicknesses ranged between 41 and 44 inches at 
locations 60, 200, and 375 yards north-northeast of the pier. Abra- 
sive action of induced currents apparently inhibited ice growth in 
zones of most vigorous flow. 


The erosive capability of water currents is manifested by recent 
experiments in the Antarctic. Specially shaped propellers driven by 
small motors were suspended through holes in the ice of McMurdo Sound. 
The propellers created vigorous currents which eroded the ice from 
below. A 10-horsepower device reportedly required 183 hours to open 
an area 30 by 85 feet in 8-foot-thick ice. An additional swath of ice 
200 feet long was eroded to a thickness of 18 inches; soon afterward, 
it fell through. 


Analysis of data obtained with the Roberts current meters revealed 
no permanent current. Mass transport of water in the area was attri- 
buted to tidal action. Peak tidal current speed was approximately 0.2 
knots (based on threshold value of the instrument). 


CONCLUSIONS 


The bubbling system operates as a huge pump capable of performing 
work on contiguous subsurface water. The rising streams of bubbles 
initiate a system of circulatory cells which extend from the bubbled 
region into adjacent water. Water from depths exceeding 15 meters is 
circulated into the agitated columns adjacent to the pier and brought 
to the surface. 


At the time of initial ice formation in 1959, convective mixing 
had occurred throughout the upper 15 meters of North Star Bugt. The 
density gradient below the 15-meter level gradually weakened with ice 
growth, and the entire water column cooled to the freezing point. 
After elimination of the warm water supply, maintenance of an open 
water area adjacent to the pier was attributed to the ice-dispersive 
and erosive activity of the induced currents coupled with the possi- 
bility that the turbulent energy also sufficed to prevent crystalline 
growth about ice nuclei. Efficiency of the system varied directly 
with turbulence. 


Considerable narrowing of the polynya by mid-winter was attrib- 
uted to marked increase of tensile strength with consequent increased 
resistance to erosive action of the induced currents as the tempera- 
ture of the ice dropped below -9.2° F. Vertical growth of the newly 
formed ice cover within the bubbled area was inhibited by this ero- 
sive action because the undersurface of the ice is weakest when its 
temperature equals that of the water. Increase in the dimensions of 
the ice-free area was observed to concur with an increase of air 
temperature to approximately O° F in early spring. This increase was 
attributed to marked decrease of tensile strength with consequent 
decreased resistance to erosion as the temperature of the ice rose 
above =-9.20 F, 


Unique properties of fresh water make the bubbling system highly 
suitable for lakes and to a somewhat lesser extent for brackish estu- 
aries. The system is less effective in salt water, because maximum 
density of water with salinity in excess of 24.7 O/oo is attained at 
the freezing point. However, factors other than the upward transport 
of warm water, as previously discussed, also contribute to the main- 
enance of an ice-free area. 


In regions where upward circulation of sensible heat is not a 
factor, maintenance of an ice-free area is predominantly dependent 
upon speed and intensity of the induced currents. 
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APPENDIX I 
MODEL OF THE POLYNYA CIRCULATORY SYSTEM 


Figure 7 illustrates streamlines generated by the motion of a 
solid sphere in an infinite mass of frictionless fluid. If we take the 
origin at the center of the sphere and the x axis in the direction of 
motion, the normal velocity (V,) at the surface of the sphere is U cos @, 
where U is the velocity of the center. 


Lamb shows that the stream function due to the sphere is 


3 
Eeibyiee, a8 
ye aU , sin 6g 


where a is the radius of the sphere, r is the radius vector from the 
center to points on or exterior to the sphere (r?a), and @ is the angle 
between the radius vector and the x axis. At any given instant the 
trajectories of the fluid particles are tangent to the streamlines. 


The total flux through a curved surface § is Wecaac 
Arbitrarily making this value equal to -2m7W, we have 


-2ny=f VndS. 


In the case where S is the surface of the above sphere (r=a) substitu- 
tion of 2@myds for dS yields 


~Qry =e Onyds, 


where ds, as shown in Figure &, is an infinitesimal length of are sub- 
tended by an infinitesimal angle, d @, on the surface §. Substitution 
of U cos @, a sin @, and ad@ for Vy, y, and ds, respectively, and 
integrating between the limits 0 and @ yields 


-y=Uo? [8cos Asin 0d 8. 
0 
Therefore, 
W=-3Ua*sin® @. 


Lamb shows that the stream function from an n pole is given by 
a-l 


Gmiicosud 
=K = 
y ex" ! 


Since the sphere acts as a dipole, 


From the boundary value r=, 


eecaee tte 2 
qin O= > Ua sin @. 
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Therefore, 


Se eee 
K= 5 Ud 


and for the general case$ i.e Ce, r>a 
3 
Sei caee 
VS> aU Si Cl. 


A model of the polynya circulatory system can be formulated from 
the idealized case by adaptation of the principles to the bubbling sys- 
tem. Considering the motion of each ascending bubble to be directed 
along the positive-downward Z axis, there will be a streamline coinci- 
dent with the Z axis and a vertical flow of water particles. Ascending 
motion, represented by a negative vertical velocity, creates divergence 
at the surface. Approximately midway between bubble streams is a 
region of convergence with consequent descending motion, clearly dis- 
cernible in Figure 9. 


Surface water beyond the pipe furthest from the pier flows out- 
ward to a distance determined by the horizontal momentum of the water 
particles. 


The data show greater density in water brought to the surface by 
the bubble activity during the pre-freezeup and initial freezeup peri- 
ods. Consequently, as the higher density surface water flowing outward 
from the divergence zone above pipe #4 suffers a gradual decrease in 
the horizontal component of the velocity vector, the vertical component 
increases. From the point where the horizontal component becomes zero, 
descending motion extends to depths where divergence directs a horizon- 
tal component toward the pier. 


The proposed model of the polynya circulatory system is presented 
in Figure 10. This cross-sectional view shows the eastern ends of the 
polyethylene pipes; arrows indicate principal paths of the water parti- 


cles. 
The author is indebted to Dr. Lloyd Simpson of the Hydrographic 


Office for advice and assistance in application of hydrodynamic prin- 
ciples in development of this idealized model of the bubbling system. 
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FIGURE 7 STREAMLINES PRODUCED BY A SOLID SPHERE 
MOVING THROUGH AN INFINITE MASS OF 
FRICTIONLESS FLUID. 


FIGURE 8 
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APPENDIX ITI 


OCEANOGRAPHIC PROGRAM - 1960 
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APPENDIX II 
OCEANOGRAPHIC PROGRAM - 1960 


A program of oceanographic data collection similar to that of 
1959 was conducted during the fall of 1960. Observations of the 
formation and growth of ice on North Star Bugt were initiated on 6 
October. On this date grease and pancake ice were observed in the 
shallow water area northeast of De Long pier. By 30 October most 
of the bay was covered with drifting floes of young ice approxi- 
mately one inch thick. Strong easterly winds with gusts to 50 
knots completely cleared the bay of ice on 3 November. Ice began 
to form again on 5 November, and a ten=tenths concentration of 
young ice was attained by 7 November with exception of an ice-free 
area adjacent to the pier. 


Commencing 8 October and terminating 7 November, serial tempera- 
ture and salinity data were obtained at 4 stations. The locations 
of Stations 1 and 2 concurred with the locations of Stations 1 and 
2 for 1959 as shown in Figure 2. Station 3 was located approximately 
100 feet north of Station 1, while Station 4 was located just off 
the shoreward end of the pier. The data are presented in Appendix 
VI. Surface temperature and salinity values for Station 1 are plot- 
ted in Figure 11. 


Data were taken at Station 4 for comparison of the oceanographic 
structure outside the bubbled area with that of the water column at 
Station 1 during the early period of ice growth on the bay. Occupa- 
tion of Station 4 necessitated breaking through the ice cover. Sharp 
rises in surface water temperature were observed on 21 and 24 October; 
easterly winds with speed maximums of 51 and 48 knots, respectively, 
were recorded on these dates, Although no data below the 10-meter 
level are available, it is evident, as indicated by the temperature 
and salinity data presented in Appendix VI, that the wind affected 
vertical mixing throughout North Star Bugt. 


On 10 October, the surface temperature at Station 1 was -1.54° c; 
the salinity was 32.30 2/00. On 15 October, the surface temperature 
at Station 2 was -1.77° C with grease ice forming in the area; sur- 
face temperature in the ice-free bubbled area was -1.689 C. Surface 
values of -1.81° C and 32.82 °/o0 were recorded at Station 1 on 29 
October; the bubbling system was not in operation, and a considerable 
amount of grease and slush ice was forming on the bay. 


Activation of the bubbling system on the following day resulted 
in quick dispersal of all ice from the bubbled area. Light grease 
and slush ice being swept from divergent regions and transported 
beneath the surface in convergent regions confirmed one aspect of 
the proposed model of induced circulation. Surface temperature of 
-1.02° C within the bubbled area indicates supercooling, since the 
calculated freezing point was -1.79° C. 
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When compared to data obtained at Station 1, those obtained at 
Station 4 on 5, 6, and 7 November indicate that vertical trans- 
port of sensible heat was not a factor in maintenance of the 
ice-free area adjacent to the pier. The temperature beneath 
the ice outside the bubbled area was identical to that of the 
isothermal water column at Station 1. 


Subsequent history of the polynya was similar to that of 
the previous winter. 
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FIGURE || SURFACE TEMPERATURE AND SALINITY AT STATION 1,8 OCTOBER- 
9 NOVEMBER I960 
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APPENDIX III 
OCEANOGRAPHIC DATA 
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APPENDIX V = SYNOPTIC METEOROLOGICAL OBSERVATIONS, 1959 
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